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ABSTRACT 


This  Research  Contribution  describes  a 
■ethodology  for  assessing  eneay  ability 
to  trail  friendly  ships  at  sea.  It  con¬ 
sists  of  four  parts.  The  first  part 
treats  the  search  for  a  lost  quarry  by 
shipborne  helicopter  or  long-range  re¬ 
connaissance  aircraft.  The  second  des¬ 
cribes  a  Markov  Model  yielding  the  frac¬ 
tion  of  ties  the  ship  Is  free  of  trill. 
The  third  part  estivates  eneay  aircraft 
requlreaents  to  achieve  specific  search 
results .  The  last  part  presents  and 
documents  an  AFL  program,  TRAIL,  that 
perfores  all  required  calculations. 
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TRAIL:  A  SHIP-TRAILING  MODEL 


This  Research  Contribution  describes  a  methodology  to  assess  enemy 
ability  to  trail  friendly  ships  at  sea.  It  consists  of  four  parts.  The 
first  part  treats  aircraft  or  helicopter  search  for  a  lost  contact.  The 
second  describes  a  Markov  model  yielding  the  fraction  of  time  the  ship 
is  free  of  trail.  The  third  part  estimates  enemy  aircraft  requirements 
to  achieve  specific  search  results.  The  last  part  presents  and 
documents  an  APL  program,  TRAIL,  that  performs  all  required 
calculations . 

SEARCH  BY  AIRBORNE  RECONNAISSANCE  PLATFORMS 

The  trailing  scenario  is  as  follows.  Each  friendly  ship  has  an 
enemy  ship  assigned  to  trail  it  wherever  it  goes.  When  trail  is  broken, 
a  helicopter  flies  from  the  trailing  ship  to  search  for  the  target.  If 
the  search  fails,  long-range  radar  reconnaissance  aircraft  are  summoned 
from  distant  bases. 

The  General  Case 


Suppose  the  friendly  ship  breaks  trail  at  time  t  -  0.  We  assume 
here  that  the  air  search  is  random  with  a  uniform  target  distribution 
(that  is,  the  target  may  be  found  anywhere  in  its  area  of  uncertainty 
with  equal  probability).  Let: 


S(t)  H  cumulative  area  searched  up  to  time  t 
A(t)  =  target's  area  of  uncertainty  at  time  t 
P(t)  =  probability  of  finding  target  before  time  t. 


Then  we  have  the  Koopman  formula: 


P(T) 


-  K  exp  /  -  ^  dt  , 


(1) 


where  K  is  determined  by  the  initial  condition.  This  formula  is 
proved  as  follows.  Let  Q(t)  Hi-  P(t).  The  probability  Q(t  +  At) 
of  not  finding  the  target  between  times  t  and  t  +  At  *  [the 
probability  Q(t)  of  not  finding  it  within  time  t]  x  [1  -  the 
probability  u  of  finding  it  between  times  t  and  t  +  At],  We  have 
u  “  [S(t  +  At)  -  S(t))/A(t)  because  of  the  uniform  target 
distribution.  These  results  lead  to  the  differential  equation 
Q'/Q  *  -S' /A,  whose  solution  is:  Q  is  proportional  to 
exp  (-/S'dt/A).  The  stated  result  follows. 
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In  particular,  let 


Tq  =  time  when  search  begins  . 


Then  S(t)  is  of  the  form: 


S(t)  -  Y(t  -  T0)S0(t)  ,  (2) 

where  Y(x)  is  the  step  function,  equal  to  1  for  x  >  0  and  equal  to  0 
otherwise.  The  Initial  condition  is: 


P(Tq)  -  min[l,S0(T0)/A(T0)l  . 


(3) 


The  derivative  of  Y(x)  is  the  delta  function  6(x),  such  that 

Jb  o(x  -  c)f(x)dx  =  f (c)  if  c  is  in  the  interval  (a,b)  and  0 
otherwise.  Applying  this  result  to  equation  1,  we  obtain: 


P(T)  -  Y(T  -  Tq)  Pq(T)  , 


where: 


PqCT)  =  1 
A  Special  Case 


Let: 


(4) 


(5) 


v  3  target  speed  after  it  breaks  trail 
v'  =  searcher's  speed 
r  =  searcher' 8  radar  detection  range. 


Assuming  that  the  target's  course  is  random,  the  target  can  be 
found  after  time  t  anywhere  inside  a  circle  of  radius  vt .  Thus: 

A(t)  -  nv2t2  .  (6) 


The  searcher  sweeps  out  equal  areas  In  equal  times,  beginning  at  TQ,  and 
so: 


S(t)  -  2rv ' ( t  -  T0)Y(t  -  Tq) 
According  to  equation  5,  we  have: 


(7) 


P0(T)  -  1  -  Clexp(c2/T) 


(8) 


with  c 


2rv' 


and  Cj 


=  exp(-c2/TQ) 


(9) 


2 

Note  that  r,  v,  and  v'  enter  only  through  the  combination  rv'/v  . 
Thus,  variations  in  r,  v,  and  v'  that  leave  rv’/v^  unchanged  do 
not  affect  the  search  results. 


A  plot  of  the  detection  probability  P(T)  versus  T  is  found  in 
figure  1.  As  can  be  seen,  P(T)  increases  very  rapidly  with  T  when 
the  search  begins  at  time  TQ.  Because  the  area  of  uncertainty  grows 


FIG.  1:  BEHAVIOR  OF  AIR  SEARCH  PROBABILITY 
DISTRIBUTION  P(T) 
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faster  than  the  searcher's  search  rate  (quadratically  versus  linearly), 
the  target  is  not  certain  to  be  found.  Indeed,  letting  T  tend  to  ® 
in  equation  8  shows  that  the  probability  of  ever  finding  the  target  is: 


Pmax  *  “  C1  ’ 


which  is  smaller  than  1, 


For  helicopter  searches  immediately  following  loss  of  trail,  we 
assume  that  the  ship's  area  of  uncertainty  expands  at  a  rate  V2  away 
from  the  trailer  and  at  a  lower  rate  v^  towards  the  trailer  (see 
figure  2).  The  difference  is  due  to  the  fact  that  the  ship  is  less 
likely  to  go  in  the  general  direction  of  the  trailer.  Assuming  the 
higher  noise  produced  by  the  accelerating  ship  does  not  reduce  its  area 
of  uncertainty,  we  see  that  we  must  simply  replace  v^  with 
(v^  +  V2)/2.  in  equation  9.  If  the  helicopter  stops  at  time  T  fron 
trail-break,  and  there  are  n0pp  independent  opportunities  per  day  to 
break  trail,  and  circumstances  propitious  for  trail-breaking  occur 
independently  every  day  with  probability  w,  then  the  daily  probability 
of  breaking  trail  is: 


-  "  [‘  -  Po 


(T  ) 
max 


X  Last  known  trailer 
position 


FIG.  2:  SHIP'S  AREA  OF  UNCERTAINTY  AT  TIME  t 


The  probability  w  can  be  tied  to  weather  conditions,  such  as  sea 
state  or  visibility,  or  to  availability  of  escorts,  or  to  other  factors  . 


Expected  Time  Target  Is  Found 


Because  the  target  is  not  certain  to  be  found,  the  expected  time 
the  target  is  found  is  infinite.  Nevertheless,  one  can  define  a 
conditional  expectation  time: 


E(T')  5  expected  time  target  is  found  given  it  is  found  before  T' . 


It  is  given  by  the  expected  time  the  target  is  found  up  to  time  T' 
divided  by  the  probability  P(T’ )  of  finding  it  before  T' ,  that  is: 


E(T' ) 


Y(T»  -  T  ) 


(T)dT  . 


The  result  is: 


E(T* ) 


-  c  exp(c„/T' 


c,  c~/T' 
-  f(— |r)  c 


k 


C2  c2^n 

+  f  (-|~)  e  L  ° 


where: 


f(c)  =  e"^ 


(l  +  £  A 

\  n»l  nnl  / 


The  value  of  f(c)  mist  be  calculated  numerically .  It  is  plotted  and 
tabulated  in  figure  3.  All  we  need  is  the  summation  term,  but  we 
introduce  f(c)  because  it  varies  much  more  gently  with  c  than  the 
summation  term.  The  latter  converges  rapidly,  especially  when  c  is 
smaller  than  about  4.  A  word  of  caution:  the  function  f(c)  is 
sometimes  needed  to  a  high  degree  of  accuracy  (several  decimal  places), 
so  that  figure  3  will  not  always  suffice.  (Example:  when  exp  (c2/TQ) 
is  very  large.)  A  TI-59  calculator  program  to  compute  f(c)  is 
provided  in  table  I.  Bablished  tables  can  also  be  used  (see  below). 


This  result  is  derived  as  follows.  The  derivative  P'  of  P  is: 


p*(T)  -  6(T  -  Tq)  Pq(T)  +  Y(T  -  Tq)  P^(T) 


-5- 


TABLE  1 


TI-59  PROGRAM  TO  CALCULATE  f(c) 


Definition:  f(c)  =  e-0 


(l  +  £  £l) 

\  n«l  nn!  / 


Listing: 


000 

42 

ST0 

021 

65 

041 

042 

08 

66 

08 

F'flU 

00 1 

10 

10 

022 

53 

< 

043 

43 

RCL 

002 

01 

•i 

1 

023 

43 

RCL 

044 

06 

06 

003 

42 

ST0 

024 

10 

10 

045 

85 

+ 

004 

06 

06 

025 

45 

t  i 
(  •• 

046 

01 

1 

005 

42 

STD 

026 

43 

RCL 

047 

95 

s 

006 

07 

07 

027 

06 

LI  t« 

048 

42 

STD 

007 

43 

RCL 

028 

54 

) 

049 

06 

06 

008 

10 

10 

029 

55 

-r 

050 

1  0 

E* 

009 

94 

+  /  - 

030 

43 

RCL 

051 

61 

GTQ 

0 1 0 

“i 

INV 

031 

06 

06 

052 

•**i 

v  r 

IFF 

011 

23 

LNX 

032 

55 

■T 

053 

76 

LBL 

012 

42 

STD 

033 

43 

RCL 

054 

10 

E* 

013 

08 

08 

034 

07 

07 

055 

36 

PGM 

014 

76 

LBL 

035 

95 

= 

056 

16 

16 

0 1 5 

O  "7 
•*»  1 

IFF 

036 

85 

+ 

057 

1 1 

R 

016 

43 

RCL 

037 

43 

RCL 

U58 

36 

PGM 

017 

10 

10 

038 

08 

08 

059 

16 

16 

0 1 8 

94 

+.-■  - 

039 

95 

= 

u60 

13 

C 

0 1 9 

INV 

040 

42 

STD 

U61 

42 

8  1  u 

020 

23 

LHH 

062 

063 

07 

92 

07 

RTN 

To  use:  Key  in  c,  and  then  press  R/S.  Display  will  flash  the  partial 
sums  for  visual  inspection  of  convergence.  9top  execution  (R/S)  when 
result  does  not  seem  to  grow  bigger.  The  result  is  in  register  8;  the 
number  of  terns  added  in  register  6.  Program  16  of  the  aaster  library 
Issued  to  calculate  n!  .  Registers  1  through  10  and  labels  A  through  E, 
E',  and  1*  are  used.  Hie  constant  c  is  in  register  10.  and  nl  is  in 
register  7.  The  larger  c  is.  the  longer  it  takes  to  get  the  result. 


(16) 


with  PQ  in  equation  8.  This  gives  (for  T'  >  TQ): 


E(T') 


P  (T  ) 
o  o 

~rrn 


i 

ttt 


/ 


TF'(T)  dT 
o 


This  is  a  general  result  for  all  P(T)  of  the  form  of  equation  4.  When 
T' ♦  T  ,  it  can  be  shown  that  E(T')  +  TQ,  as  expected.  It  can  also  be 
shown  that  E(T')  is  between  TQ  and  T' ,  as  expected.  In  this 
special  case,  the  first  term  is  0  because  P0(T0)  ”  0  (from  equation 
8).  Expanding  the  exponential  in  the  integrand  term  by  term  gives: 


TP'(T)  dT 
o 


(17) 


and  the  stated  result  follows.  The  sum  can  be  related  to  the 
exponential-integral  function  Ei(x)  by  [1,  p.  229,  5.1.10]: 


«  n 

22  — f  ”  Ki(c)  -  Y  -  In  c  ,  (18) 

n*l 

where  Y  2  .5772  is  Catalan's  constant  and  Ei(c)  is  tabulated  [1,  pp .  238 
ff]. 

When  P(«)  is  close  to  1  (that  is,  when  the  searcher  is 
practically  certain  to  find  its  target)  it  is  tempting  to  let  T'  tend 
to  <»  in  equation  13  to  get  Che  expected  time  the  target  is  found. 
Unfortunately,  E(T' )  diverges  because  of  the  In  (T' )  term.  However, 
it  is  a  mild  divergence,  so  that  letting  T'  be  a  large  but  reasonable 
constant  gives  the  desired  expected  time,  and  the  latter  does  not 
Increase  ouch  with  large  increases  in  T' . 


Resumption  of  trail 

Once  the  aircraft  finds  the  ship,  only  half  the  work  is  done.  The 
trailing  ship  oust  now  resume  trail.  What  was  the  trailer  doing  after 
trail  was  broken?  Because,  in  the  long  run,  any  course  the  trailer 
takes  is  equally  likely  to  be  toward  the  final  target  position  (at  the 
time  the  aircraft  finds  it)  as  away  from  it,  we  assume  the  trailer 
roughly  maintains  its  initial  position  0  (see  figure  4) .  At  the 
time  T  the  ship  is  found,  the  area  of  uncertainty  has  grown  to  a 
circle  (C)  of  radius  vT  centered  at  0.  The  ship  is  found,  on 
average,  at  a  point  a  distance  R_  from  0  such  that  the  circle 
(CQ)  of  radius  RQ  contains  half  the  area  of  uncertainty.  This  gives: 


Rq  «vT//r  . 


(19) 
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Ship  area  of 
uncertainty 
(C) 


Point  where  aircraft 
finds  ship 


FIG.  4:  TYPICAL  GEOMETRY  WHEN  AIRCRAFT  FINDS  SHIP 


Once  Che  aircraft  communicates  the  ship's  position  to  the  trailer, 
the  trailer  moves  toward  the  ship.  But  the  ship  moves  too.  In  the  long 
interval  after  trail  is  broken,  the  ship  might  receive  intelligence 
information  on  the  trailer's  movements,  or  at  least  an  estimate  of  the 
trailer' 8  position  when  the  aircraft  detected  the  ship,  and  it  could 
move  away  from  that  position.  If  so,  it  may  be  a  long  time  before  the 
ship  is  back  under  trail,  if  it  ever  is.  Assume  a  ship-quarry  closure 
rate  of  Vj.  Then,  if  the  aircraft  detects  the  ship  at  time  T,  the 
ship  is  back  under  trail  at  time  kT,  where: 


k  = 


1  +  vt/VT 

v*. 


(20) 


More  than  one  aircraft  may  be  needed  for  any  one  mission,  with  each 
aircraft  relaying  its  findings  to  its  relief .  Note  that  the  time 
between  T  and  kT,  when  the  aircraft  is  holding  the  ship  and  waiting 
for  the  trailer  to  arrive,  can  in  some  sense  be  considered  a  time  under 
trail  for  purposes  of  possible  attack  by  platforms  other  than  the 
trailing  ship.  In  our  context,  however,  "trail"  refers  only  to  direct 
trail  by  a  specific  trailing  ship  unassisted  by  outside  assets. 


Expected  Kate  of  Ifrail  Reacquisition 

Our  model  assumes  a  constant  probability  per  unit  time  that  a 
trailer  will  reacquire  Its  target  ship  after  trail  Is  broken.  Urns, 

If  s  Is  the  probability  of  trail  being  resumed  nj  days  after  It  Is 
lost,  then  the  dally  trall-reacqulsitlon  probability  P'  Is  such  that: 

1  -  (1  -  P')"1  -  s  ,  (21) 


the  days  being  assumed  Independent.  Solving  the  equation  for  P' 
when  s  Is  known,  we  get: 

1/n. 

P'  -  1  -  (1  -  s)  .  (22) 


here,  the  ship  Is  first  reacquired  by  an  aircraft  at  time  T, 
distributed  according  to  equation  4.  After  the  aircraft  has  found  the 
ship.  It  stays  with  it  until  the  trailer  can  arrive  and  resume  trail, 
which,  as  described  above,  takes  a  period  of  time  proportional  to  T. 
Thus,  let: 


k  =  constant  such  that  trail  Is  resumed  at  time  kT  if 
aircraft  finds  ship  at  time  T. 

sQ  =  ratio  of  desired  unit  of  time  for  reacquisition  rate  P' 
to  unit  of  time  in  distribution  P(T)  (here,  sQ  ■  24 
since  we  want  a  rate  per  day  and  T  is  in  hours)  . 


q  2  sQ/k  . 


What  we  want  is  the  expectation  value  P'  of: 


1  -  [1  -  P(T)]<l/T 


(23) 


relative  to  the  distribution  P(T)  in  equation  4.  For  the  reasons 
explained  earlier,  we  assumed  that  the  search  will  be  called  off  at  a 
certain  elapsed  time  t  .  This  time  is  designed  to  ensure  that  the  ship 
will  be  found  with  a  given  probability  p,  smaller  than  1  -  cj  .  t  is 
therefore  given  by  inverting  equation  8: 
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Ue  have: 


?’  “  £  (  1  -  f1  -  p  (T)]q/T  (  P'(T)  dT 


Using  equation  15,  we  have: 


P’  -  1 


p.d.) 


,/T0  Po(To) 
PC(T) 


/  [l  -  po(t>]',/t  p;  <t>  iT  • 


which  Is  valid  Cor  all  PQ .  Specializing  to  the  distribution  in 
equation  8,  we  get: 

P.  .  1  _  Clq/To  exp(c2q/T^  [l  -  c^xp^/T^J  P^Ct) 

-p"1  (t)  £  Cq/T  |exp(c2q/T2^J  c^T**2  exp(c2/T)dT  .  (27 

o 

Specializing  further  using  equation  9,  which  indicates  that  PQ(T0)  -  0 
we  have,  with  the  change  of  variable  x  ■  1/T, 

P'  -  1  -  l'\*  exp  [!£  x2  ♦  (c2+  !g  1„  .,)  ,]  .  m 

This  integral  must  be  calculated  numerically  (see  table  2  for  an  APL 
program).  Note  that,  unlike  the  expected  time  to  find  the  target  we 
calculated  earlier,  P'  does  tend  to  a  limit  when  t  ♦  •». 

FRACTION  OF  TIME  FREE  OF  TRAIL 

The  primary  measure  of  effectiveness  in  avoiding  trail  is  the 
fraction  of  time  at  sea  the  ship  is  free  of  trail. 

Result 

We  now  describe  the  Markov  model  that  leads  to  the  following 
result: 

E  =  fraction  of  time  at  sea  the  ship  is  free  of  trail 

p  .  0"  -  V'  -  pop)  (1  -  p  -  P')  [i  -  (1  - p  -  P')"] 

“  p  +-51  + - 5 -  * 

P  n(P  +  P'/ 


TABLE  2 


NUMERICAL  INTEGRATION® 


Listing: 

VXMTEGRALCm? 

7  INTfAS  INTEGRAL  ERSJHJHf K)L 

C13  w«-if0,5x(  (k*-0)++/sntegranp  a»)*H(.- /#a» 

£23  r  :*»«•«»  <**c<f«) -X- 134-2)  +«*+/! ntegrand<  ifA»)  +  (H«-H*2 >x “1+2X 
\2*'’l+K«-|<  +  l+,-<-0 

£33  o :»«.»»,  <  <4tL)xwEfW3>-wC<f**>-K3>*"l+4*i-*.L+i 

C43  40xiL<k 

C53  -»rxiERS<  I -/**£ (f **«■  (-l+2xX)t»)-0»K3 

V 

b 

How  To  Run:  To  calculate  J  f(x)  dx,  define  the  function  f(x) 

a 

in  APL  (calling  it  INTEGRAND) ,  input  the  two-vector  AB  - 
(a,b)  and  the  desired  precision  EPS,  and  then  key  in  AB 
INTEGRAL  EPS.  Vector  origin  must  be  I  (□  10  ♦  1). 

Example :  a  -  0,  b  ■  1/8.4,  EPS  *  10-^, 

f(x)  -  exp  [239. 5x^  +  (38.2  +  6.27  In  Cj)  x] 

V  INTEGRAND  [[]]  7 
v  t*- integrand  >: 

113  rt-A  (23?.5xx*2)+Xx33,2+6.27xaci 
v 

□  104-1 
cl4-.010d 

(  0  f  -r-8  ♦  4  )  INTEGRAL  1E~6 

1 ♦ 6409036 


a  Using  Romberg's  rule.  Program  is  from  CNA's  APL  Library. 
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PQ  =  probability  the  ship  is  trailed  when  it  leaves  port 


P  =  probability  the  ship  will  shake  free  on  a  day  when 
it  starts  out  trailed  (assumed  constant) 

P'  =  probability  the  ship  is  reacquired  by  its  assigned 
trailer  on  a  day  when  it  starts  out  free  (assumed 
constant) 

n  =  number  of  days  the  trip  lasts. 


The  first  term  frequently  dominates.  A  required  correction  to  E 
will  be  discussed  shortly. 

Assumptions 

The  time  unit  is  arbitrary.  We  chose  to  use  the  day.  This  means 
that  in  a  given  day,  the  ship  is  assumed  to  be  either  trailed  all  the 
time  or  free  all  the  time.  In  other  words,  we  assume  that  the  process 
of  breaking  trail  or  reacquiring  trail  occurs  instantaneously  at  the 
beginning  of  a  day.  This  simplifying  artifice  does  not  affect  the 
result.  If  the  unit  of  time  is  chosen  to  be,  say,  the  hour,  then  the 
"rates”*  P  and  P'  must  be  divided  by  24**  and  the  time  interval  n 
multiplied  by  24.  It  can  be  verified  that  the  resulting  E  remains 
virtually  unchanged. 

Note  that  the  dally  break-trail  probability  P  refers  to  losing  trail 
for  a  large  fraction  of  a  day.  Broken  trails  that  are  quickly 
reacquired,  say  by  helicopters  operating  from  trailing  ships,  are  not 
counted.  Equation  11  gives  a  possible  estimate  of  P. 


*  Strictly  speaking,  P  and  P'  are  transition  probabilities,  not 
rates — a  subtle  but  important  difference  in  statistical  calculations. 

If  P  were  a  rate,  there  would  be  a  probability  e  -pPn/n!  of  n 
transitions  per  day;  in  fact,  there  can  be  only  one  transition  per  day, 
and  its  probability  is  P,  not  Pe“.  However,  since  the  probability 
of  a  transition  in  n  days  is  1  -  (1  -  P)n,  and  since  this  expression 
can  be  rewritten  1  -  exp  (-Xn)  if  one  wishes  to  consider  it  a  Poisson 
process,  one  can  identify  -ln(l  -  P)  with  the  "rate"  X.  For  small 
P,  X*P.  . 

**  Or  replaced  by  P  -  1  -  (1  -  ?)1'  ,  assuming  independence  from  hour 

to  hour. 
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It  is  also  assumed  that  each  day  is  independent  of  the  previous  day  for 
breaking  and  reacquiring  trail. 

Range  of  Validity  of  Model 


Obviously,  the  closer  the  exact  probability  of  reacquiring  (or 
breaking)  trail  over  A  days  resembles  f(A)  ■  1  -  (I  -  p)A,  where  p 
is  a  constant,  the  better  the  model.  Note  that  the  probability  of 
reacquisition  is  0  for  some  time  and  then  surges  to  high  values  in  a 
short  period  of  time,  unlike  the  even  behavior  of  f(A>*  (see  figure 
5).  Nevertheless,  the  model  is  still  valid,  provided  both  the  total 
trip  duration  n  and  the  mean  time  between  trail-breaking  events  are 
larger  than  the  time  Aq  when  the  two  distributions  do  not  differ 
significantly  in  value.  If  these  conditions  are  not  met — say  if  trail 
is  frequently  broken — then  the  frequent  Initial  periods  of  no-trail  will 
add  up  significantly,  and  the  model  will  underestimate  the  fraction  of 
time  the  ship  is  free  of  trail. 


FIG.  5:  ILLUSTRATION  OF  APPROXIMATION 
IN  THE  TRAIL  MODEL 


*The  exact  probability  of  trail  resumption  within  A  days  after  trail 
is  broken  is  0  if  A<kT0/s0  and  P0(s0A/k)  otherwise,  with  PQ  as  in 
equation  8. 
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The  Markov  Chain 


Since  the  days  are  assumed  independent,  we  have  a  Markov  process, 
illustrated  in  figure  6.  It  has  two  states:  free  (F)  and  trailed  (T)  . 
Ihe  state  vector  is  S  -  (Pp.Pj),  where  PA  is  the  probability  of  being 
in  state  A.  Ihe  initial  state  SQ  is: 

SQ  -  (l  -  P0*po)  *  OO) 


ETC. 


FIG.  6:  MARKOV  CHAIN  FOR  TRAILING  MODEL 

If  P^g  is  the  daily  probability  of  transition  from  state  A  to  state  B, 
then: 


PpF  m  1  “  p' »  **FT  “  p'  »  PTF  ”  p*  PTT  "  ^  ”  P 


(31) 


The  transition  matrix  M  is  therefore: 


(32) 
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JtVi .  - 


The  state  vector  after  k  days  is: 


\  “  soM  • 


(33) 


The  kth  power  of  the  matrix  M  can  be  calculated*  to  be: 


1  /  P  +  P'W  P*  -  P'w\ 

+  P'  \  P  -  PW  P'  +  PW  )' 


(34) 


where: 


W  =  (1  -  P  -  P’ )K  • 


(35) 


*To  find  the  nth  power  of  a  2x2  matrix  A,  assume  An  ■  m^A  +  m2l, 
where  I  is  the  identity  matrix.  If  A  ^  and  A2  are  the  eigenvalues 
and  X  an  eigenvector  of  A,  we  have  AX  *  XX  and  AnX  -  AnX,  which 
gives  An  -  mjA  +  m2  for  A  =  Aj  or  A  -  A2,  and  yields  mL  and  m2 
in  terms  of  A^  and  A2.  The  result  is: 

ml  “  cn  and  m2  "  -xlx2  Cn-1  • 


which  gives: 


/  aCn  "  XlX2Cn-l 

V  CCn 


dC 

n 


\ 

"  XlX2Cn-l/ 


> 


where: 


A 


The  eigenvalues  are  obtained  by  writing  that  the  determinant  of  A  -  A I 
is  0.  This  yields  A^  ■  p  +  q  and  A2  =  p  -  q,  where: 

p  -  5-+-!  ,  q  -  [(d  _  a)2  +  4bc]l/2  . 


In  our  case,  Aj  ■  1  and  A2  -  1  -  P  -  p' . 
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The  probability  that  the  ship  is  free  on  the  kth  day  is  the  first 

element  of  S^,  namely: 


-  (P  +  P')-1  [P  +  W(P'  -  PQP'  -  P0P)]  . 


(36) 


The  expected  nunber  of  days  free  of  trail  is  simply  the  sum  of  the 
probabilities  of  being  free  for  each  of  the  days  at  sea.  Dividing  by 
n  gives  the  untrailed  fraction  E: 


E 


n 


n 


L 

k-l 


(37) 


which  leads  directly  to  the  stated  result  after  use  of  the  geometrical 
progression: 


n 


E 

k-0 


1  -  x"+1 

1  -  X  * 


(38) 


Special  Cases 

When  P'  *  0,  which  is  generally  the  case  when  the  trailer  is  left 
to  its  own  devices,  the  fraction  E  of  time  free  of  trail  simplifies  to 
E  -  EQ,  where: 


E 

o 


1  - 


?G  (1  -  P) 
nP 


(39) 


*  This  special  case  can  be  simply  derived  as  follows.  The  probability 
of  shaking  the  trailer  off  on  the  kth  day  is  P0(l-P)k_*P.  This  action 
results  in  n  -  k  +  1  days  free  of  trail.  Thus: 

n  k-l 

nE  -  n(  1  -  P  )  +Z  (n  -  k  +  1)P  (1  -  P)*  P. 
o  o  .  ,  o 

k-l 

Using  the  geometrical  progression  result  of  equation  38  and  the  fact 
that 

E  kxk-  x  4—  E  ■  x(l  ~  x)  ^  [l  -  xn  -  n(l  -  x)xnl  , 
k-0  dx  k-0  L  J 

we  get  the  stated  result . 
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When  n  -*•  ®  (long  trips),  we  have: 


E 


P 

+  P  +  P' 


(AO) 


When  P  =  P'  (the  ship  breaks  trail  as  easily  as  the  trailer 
reacquires)  and  P  -  V2  (the  ship  is  trailed  during  port  exits  as  often 
as  not),  the  ship  Is  trailed  exactly  half  the  time.  The  same  result 
obtains  if  P  *  P'  *  V2  ,  or  if  P  ■  P'  and  n-*»,  or  if  P  ■  P*  *  PQ  ■  1 
and  n  is  even. 

Behavior  of 

From  equation  36  we  infer  that  Q^,  the  probability  of  being  free 
on  day  k,  decreases  with  k  if  P'  -  PQP’  -  PQP  is  positive  and 
increases  with  k  if  that  quantity  is  negative,  provided  P  +  P'  is 
smaller  than  1.  The  asymptotic  limit  (steady-state  value)  of  for 

large  k  is  P/(P  +  P').  If  P  +  P'  is  larger  than  1,  oscillates 

up  and  down  around  the  limit,  but  the  oscillations  damp  out  with  time. 

If  P  +  P'  *  1,  Qk  is  always  equal  to  the  limit. 

Examples 

The  behavior  of  is  illustrated  in  figure  7  for  P  **  P* , 

PQ  *  .2,  and  n  ■*  90.  As  can  be  seen,  the  closer  P  +  P'  is  to  ., 
the  faster  the  process  reaches  its  steady  state.  The  latter  is  reached 
rather  fast  anyway:  by  the  tenth  day,  differs  from  the  steady- 

state  value  by  less  than  10  percent. 


FIG.  7:  DAILY  PROBABILITY  Q*  SHIP  IS  FREE  OF  TRAIL 
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Figure  8  shows  how  the  fraction  of  time  free  of  trail,  E,  varies 
with  P,  P',  and  PQ.  For  simplicity  and  so  as  not  to  introduce  bias, 
we  choose  P  -  P' — the  ship  breaks  trail  as  easily  as  the  trailer 
reacquires  it.  The  trip  lasts  n  *  90  days.  The  dependence  on  PQ, 
predictable  in  nature  but  surprising  in  magnitude,  is  as  follows:* 

•  If  P  is  relatively  large  (over  5  percent  a  day) — 
frequent  trail  loss  and  reacquisition — the  result  is 
virtually  independent  of  PQ.  For  example,  for  P  -  5 
percent,  if  the  ship  is  picked  up  in  port  (PQ  ■  1),  it  is 
free  of  trail  45  percent  of  the  time;  if  the  ship  is  not 
picked  up  in  port,  it  is  free  of  trail  55  percent  of  the 
time. 

•  If  P  is  relatively  small  (less  than  5  percent  a  day) — 
rare  trail  loss  and  reacquisition — the  result  depends 
strongly  on  PQ.  For  example,  for  P  =  1  percent,  if  the 
ship  is  picked  up  in  port  (PQ  ■  1),  it  is  free  of  trail  27 
percent  of  the  time;  if  the  ship  is  not  picked  up  in  port, 
it  is  free  of  trail  73  percent  of  the  time. 

Corrected  Fraction  of  Time  Free  of  Trail 

If  the  aircraft  search  succeeds  only  with  probability  p,  then  equation 
28  gives  the  corresponding  P*  and  equation  29  the  corresponding  E. 
Since  the  ship  is  not  reacquired  with  probability  1,  and  since  only  98 
percent  of  the  trail-breaking  events  are  prosecuted,  the  (corrected) 

expected  fraction  of  time  free  of  trail,  E,  is: 


E  -  gp  x  E  (calculated  with  the  P'  corresponding  to  p, 

using  equation  28)  (40a) 

+  (1  -  gp)  x  E  (calculated  with  P'  ■  0,  from 
equation  39) 


where  g  -  .98. 


*  The  curves  are  symmetric  with  respect  to  the  line  E  ■  V2  •  That  is, 
if  P  ■  P*,  E  -  V2  changes  sign  when  PQ  changes  into  1  -  P  . 


-19- 


FIG.  8:  EFFECT  OF  TRAILING  AT  PORT  EGRESS 
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Equation  40a,  which  is  an  approximation,  is  a  correction  to  E 
required  by  the  observation  that  a  constant  daily  reacquisltlon 


probability 


however  small,  implies  eventual  reacquisition,  whereas 


in  reality  reacquisition  is  not  certain:  its  probability  can  be  no 
higher  than  gp.  Equation  40a  is  a  slight  underestimate  of  the  fraction 
of  time  free  of  trail. 


ENEMY  AIRCRAFT  REQUIREMENTS 

In  this  section,  the  number  of  reconnaissance  aircraft  the  enemy 
requires  to  achieve  specific  search  results  is  calculated. 

To  estimate  the  required  aircraft  inventory  for  the  reacquisltlon 
efforts,  the  quantity 


fffl(k)  =  probability  of  m  trail-breaking  events  in  the  first 
k  days  of  the  n-day  trip 


is  needed.  It  is  not  easy  to  calculate. *  Here  we  only  present  an  exact 
line  of  approach  and  an  approximation,  and  calculate  the  mean  of  the 
distribution. 

Exact  Approach 


This  method,  developed  by  James  K.  Tyson,  consists  of  defining  the 
matrix: 


\ s  (l »!',). 


which  reduces  to  the  Markov  matrix  M  in  equation  32  for  z  ■  1.  By 
"flagging"  the  P  with  a  z,  we  isolate  a  trail-breaking  event.  The 
coefficient  of  zm  in  any  result  based  on  Mz  is  then  associated 
with  m  trail-breaking  events.  Let  =  (F^  Tk)  ■  SQM^  be  the  state 
vector  after  k  days.  Then  fm(k)  is  simply  the  coefficient  of  zm  in 


*  fm(k!)  is  not,  as  one  might  think,  the  binomial  distribution 

(m)  _p)h~m.  jt  depends  on  PQ  and  P*  as  well.  As  an 

Illustration,  if  PQ  ■  P'  -  0,  the  ship  is  never  trailed,  so  it  never 
breaks  trail,  and  f_(k)  -  0  for  all  nonzero  m  and  k. 


Fk  +  Tk  (we  8um  the  two  because  it  does  not  natter  whether  the  end  state 
is  free  or  trailed),  that  is: 


f  00 

ni 


(K  + 


V  , 


-o* 


The  eigenvalues  of  Mz  are: 


(42) 


Xj  ■  v  +  i/5T  and  ^  ■  v  -  /T 


(43) 


where: 


v  =  1 


P  +  P' 
2 


We  have: 


(44) 


«,>“  ■  V1.  -  hVW  •  <«> 

where: 


where  5^  indicates  that  the  summation  extends  only  to  odd  S.  Note  that 
fm(k)  -  0  for  m>k/2.  The  formula  checks  out  for  small  values  of  k, 
but  no  simple  closed  form  could  oe  found. 

Expected  Number  of  Trail-Breaking  Events 

The  mean  of  the  unknown  distribution  fa(k)  can  be  calculated 
exactly  using  the  results  from  the  Markov  approach.  We  have: 


F(k)  =  expected  number  of  trail -breaking  events  in  the  first  k 
days 

-  P  [n(  1  -  E)  +  PQ  -  (1  -  Q^]  ,  (47) 


with  E  as  in  equation  29  and  as  in  equation  36.  The  first  term 

generally  dominates,  especially  for  large  n  and  small  E. 
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This  result  is  derived  as  follows.  The  probability  of  a  trail¬ 
breaking  event  on  day  1  is  PQP;  on  day  2  it  is  (1  -  Qj)P;  and  on  day 
k  It  is  (1  -  Qk-i)p*  (The  probability  of  being  trailed  on  day  i  is 
1  -  Q^).  Adding  up  these  probabilities  gives  F(k).  Using  the 
definition  of  E  in  equation  37,  we  get  the  stated  result. 

Approximation 

He  base  this  approximation  on  equilibrium  (steady  state)  being 
reached  rapidly.  According  to  equation  36,  the  probability  1  -  of 

being  trailed  on  day  k  rapidly  approaches  its  steady-state  value, 


1  "  \  *  P  +  P'  * 


(48) 


when  any  of  the  following  three  conditions  is  met: 


•  PQ  *  P'/(P  +  P')  (the  trip  happens  to  (49) 

begin  in  the  steady  state) 

•  P  +  P  »  1  (the  system  falls  in  the  steady  (50) 

state,  *  P,  almost  from  day  1) 

•  (1  -  P  -  P')k  «  0  for  k  «  n  (fast  approach  (51) 

to  the  steady  state). 

If  one  of  these  conditions  obtains,  we  can  assume  that  the 
probability  s  of  breaking  trail  on  any  given  day  in  the  steady  state. 


P'P 
P  +  ¥* 


(52) 


is  independent  from  day  to  day.  The  probability  f_(n)  m  trail- 

breaking  events  in  the  n  days  is  therefore  approximated  by  a  binomial 
distribution: 

fm(n)  >(”)  sm(1  -  s)n^  .  (53) 

Now  compare  the  mean  of  the  above  distribution,  ns,  with  the  exact 
mean  F(n)  calculated  earlier  (equation  47).  If  the  two  are  in 
agreement,  the  approximation  is  probably  good.  Generally,  agreement 
will  be  observed  if  both  these  two  conditions  are  met: 


-23- 


•  E  *  P/(P  +  P') 


•  n  »  PQ  (1  +  P/P')  -  1  . 


The  first  condition  obtains  for  large  n  or  if  either  of  the  condi¬ 
tions  in  equations  49  and  50  is  met.  The  second  expresses  the  fact  that 
the  first  term  should  dominate  in  the  expression  of  I(n)  in 
equation  47 . 

Mote  that  the  steady-state  distribution  of  trail  reacquisition 
events  is  the  same  as  that  for  trail -breaking  events,  as  it  should  be. 
Indeed,  PP'/(P  +  P' )  can  be  interpreted  as  either  Px[P’/(P  +  P' )]  or 
P'x[P/(P  +  P')J . 

An  APL  program  for  calculating  the  cumulative  binomial  is  given  in 
table  3. 

Calculation  of  Aircraft  Requirements 

Given  the  validity  of  this  approximation,  the  standard  deviation 
o  is: 


o  *  [ns( 1  -  s) ] 


The  cumulative  distribution,  or  probability  that  k  or  fewer  events 
will  occur,  is  approximated  [2]  by: 


1/2  -  ns 


where: 


<fr(x)  =  /  (2ir )  1/,2exp(-t2/2)dt 


is  the  standardized  normal  distribution.  Ihus,  the  probability  of  No 
events  over  the  mean  ns  is  given  by: 


For  2o  confidence  (N  -  2),  the  probability  is  over  98  percent,  since 
<p(2)  -  97.7.  Ihus,  by  planning  for  2o  events  over  the  mean,  one  is 
confident  of  being  able  to  handle  at  least  97.7  percent  of  the  events. 
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TABLE  3 


CUMULATIVE  BINOMIAL  APL  PROGRAM 


Coding: 

9»«N£D3* 

<9  >IN 

£13  *«•(****>**♦** 

£23 

£33  *r  m  *}R|'  »  RR  »  •  fRRj'  t  so  s  n  *js 

£43  -M  *  ■  |H 

£53  aj*i«.w,v*.v+<m;m)  x(S*h>  *<i-s)  tN-M«.M+i 

£63  4(M<m)/a 

£73  'FKOiABIUITT  Of  K  EVENTS  OR  LESS  (K  r»  0  r  1  *  2 »  .  .  .  t H )  1  MU" 
£83  a  MEAN  NUMBER  OF  EVENTS }  '}NX* 

£93  'STANDARD  DEVIATION  }  •  } ( N X * *  1 -S ) *0 . 5 

V 

Input :  N  (number  of  days  n),  P  (probability  P),  PP 
(probability  P') 


Intermediate  output:  the  binomial  probability  S  s  PP'/(P  +  P') 
To  execute:  Key  in  Inputs,  then  BIN. 


Example : 


*♦■♦243 

***.28 

N*90 

BIN 


key  in 


*  -  0.243  »  **  =  0.28  »  so  s  =  0.1300956 
n  «  90 


PROBABILITY  OF  K  EVENTS  OR  LESS  (K  s  0 » 1 f 2 f . . . » N > } 

3.568103E-6  0.0000515935  0.0003712048  0.001773291 

0.006333921  0.01806518  0.04291957  0.08752375 


0.1567315 

0.2510327 

0.3652662 

0.489512  0.6118379 

0.7216022 

0.811887 

0.8802984 

0.9282563  0.9594764 

0.9784119 

0.989143 

0.9948403 

0.9976804  0.9990125 

0,9996015 

0 ♦ 9998474 

0.9999445 

0.9999808  0.9999937 

0.999998 

0.9999994 

0.9999998 

11111111 

1111 

1  1  1 

1111 

11111111 

1111 

1  1  1 

1111 

11111111 

1111 

1  1  1 

1111 

1  1  1 

MEAN  NUMBER  OF  EVENTS*  11.7086 
STANDARD  DEVIATION}  3,191452 
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The  number  of  trail-breaking  eventa  per  n-day  trip  that  are  observed  98 
percent  of  the  time  is  therefore: 


f(n)  +  2/ns(l  -  s) 


(60) 


for  a  single  ship.  For  n'  independent  ships,  the  mean  is  multiplied 
by  n'  and  the  standard  deviation  by  /»',  *  so  that  the  number  is: 


v  *  n'f(n)  +  2  /nn ' s ( 1  -  8 ) 


(61) 


It  is  always  best  to  use  aircraft  assets  to  prosecute  many  trail- 
breaking  events  with  few  sorties  per  event  rather  than  few  trail- 
breaking  events  with  many  sorties  per  event,  because  the  area  of  uncer¬ 
tainty  grows  faster  than  the  area  searched. 

A  conservative  estimate  (to  be  discusi ed)  of  the  number  N„  of 
reconnaissance  aircraft  required  in  Inventory  to  guarantee  (1)  that  over 
98  percent  of  events  will  be  prosecuted  and  (2)  that  trail  will  be 
resumed  after  each  event  with  probability  p,  is: 


NR  *  Ws  •  (62) 

where: 

ns  -  number  of  sorties  an  aircraft  makes  per  month 

t&  -  maximum  time  an  aircraft  can  stay  on  station 

tj  *  aircraft  time  on  station  that  guarantees  trail  is  resumed 
with  probability  p  (»  kt  -  TQ) 

30  -  - 

v  -  —  v,  with  v  in  equation  61  (63) 

-  number  of  trail-breaking  events  per  month  aircraft  inventory 
is  tailored  to  prosecute. 


*  When  n  independent  random  variables  are  summed,  adding  the  means  and 
the  variances  gives  the  mean  and  variance  of  the  sum. 
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This  number  is  derived  as  follows.  Since  ti  aircraft -hours  are 
required  per  event  and  our  aircraft  can  only  spend  tg  hours  on  sta¬ 
tion,  n*  ■  t./tfl  sorties  are  needed  per  event,  that  is,  vnA  sorties 
per  month.  Since  each  aircraft  can  make  only  ng  sorties  per  month, 
vt»A/ns  aircraft  are  needed. 

If  each  reconnaissance  aircraft  requires  q  tankers,  the  total 
number  of  aircraft  required  is: 

Nrt  -  (1  +  q)\  .  (64) 

The  result  in  equation  62  for  enemy  aircraft  requirements  is  only  a 
rough,  conservative  estimate.  Many  more  assumptions,  models,  tradeoff 
analyses,  and  operational  data  are  required  for  an  accurate  assessment 
of  those  requirements.  First,  the  trail-breaking  events  are  assumed 
evenly  spread  in  time,  even  though  they  may  tend  to  cluster.  For 
example,  if  trail  tends  to  be  broken  in  bad  weather,  then  bad-weather 
periods  will  contain  many  more  events  than  good -weather  periods.  These 
surge  requirements  will  drive  up  enemy  requirements.  Consider  the 
example  (arbitrary  parameters)  where  five  sorties  are  needed  per  event 
for  a  p  ■  .6  reacquisition  probability  and  the  enemy  plans  for 
v  -  85  trail-breaking  events  per  month  under  conditions  yielding  a 
w  -  .12  dally  probability  that  circumstances  will  be  propitious  for 
trail -breaking.  The  aircraft  fly  12  sorties  per  month.  Our  results 
show  that  the  enemy  needs  85  x  5/12  -  35  aircraft.  But  if  we  assume 
perfect  correlation  between  the  times  when  circumstances  propitious  for 
trail-breaking  occur  throughout  the  ships'  area  of  deployment,  trail- 
breaking  events  will  cluster  every  1/.12  *  8.3  days.  Thus,  as  many  as 
85  x  8.3/30  -  24  events  could  be  expected  to  occur  roughly  at  the  same 
time,  which  would  drive  enemy  requirements  up  to  24  x  5  -  120 
aircraft — over  three  times  the  number  we  calculated.  Therefore,  the 
detailed  distribution  of  trail-breaking  events  in  time  and  Information 
on  the  correlation  between  the  occurrence  of  circumstances  propitious 
for  trail-breaking  in  different  areas  are  needed.  Second,  we  assumed 
the  enemy  knows  when  to  expect  reacquisition,  when  in  reality  he  does 
not  and  must  automatically  send  follow-up  sorties  to  ensure  that  an 
aircraft  about  to  leave  station  is  relieved.  On  the  other  hand,  we  used 
the  maximum  (not  average)  number  of  sorties  consistent  with  a  given 
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probability  of  reacquisition  (p),  because  there  is  a  limit  to  the 
extent  to  which  unused  sorties  can  be  “saved"  for  future  use.* 


Information  on  how  frequently  an  aircraft  can  be  flown  in  a  given  time 
interval,  however  short,  must  be  included  in  more  refined  esimates  of 
enemy  requirements.  Also,  variations  that  might  be  attractive  for  the 
enemy — such  as  flying  more  than  one  sortie  soon  after  trail  is  broken, 
while  the  area  of  uncertainty  is  still  small — must  be  assessed. 

As  can  be  expected  the  dominant  factor  by  far  in  anticipating  enemy 
aircraft  requirements  is  w  —  the  fraction  of  time  opportunities  for 
trail-breaking  occur. 

APL  PROGRAM  FOR  TRAIL 

Description 

Most  of  the  above  calculations  were  encoded  in  an  APL  program.  The 
26  input  variables  are  listed  in  table  A  and  the  23  output  variables  in 
table  5.  The  coding  is  in  table  6.  There  are  four  different  functions: 

•  ARM,  the  overall  control  function  (for  jiircraft 
reconnaissance  mission) 

•  INTEGRAL,  the  function  that  performs  the  integral  in 
equation  28 

•  INTEGRAND,  the  integrand  of  equation  28 

•  PEE,  the  function  that  calculates  P  from  equation  11. 

The  program  is  executed  by  keying  in  all  input  variables  and  then  ARM. 
Should  the  input  probability  of  finding  the  ship,  p,  be  larger  than 
can  possibly  occur,  an  "invalid"  message  is  printed.  Vector  origin  is 
1,  and  results  are  printed  to  four  significant  figures,  as  instructed  in 
line  2  of  ARM.  Should  the  input  tolerance  EPS  used  to  calculate  the 


*  The  following  simple  example  shows  that  there  are  cases  when  only  the 
maximum  number  of  aircraft  will  satisfy  the  requirement.  Assume  there 
is  exactly  one  trail-breaking  event  per  day,  that  the  aircraft  can  each 
fly  exactly  one  sortie  per  day,  and  that  the  ship  is  reacquired  the  same 
day  it  breaks  trail,  if  at  all.  Assume  further  that  one  sortie  yields  a 
.A  probability  of  reacquisition,  two  sorties  .6,  and  three  sorties  .8. 
The  average  number  of  sorties  flown  is  then  (.A  x  1)  +  (.2  x  2)  + 

(.A  x  3)  **  2.  The  question  is:  how  many  aircraft  in  inventory  will 
guarantee  that  80  percent  of  the  trail-breaking  events  are  reacquired 
over  a  long  period  of  time?  The  answer  is  three  (the  maximum),  not  two 
(the  average).  If  only  two  aircraft  are  bought,  there  will  be  no  way  a 
third  can  be  found  when  needed.  Unused  sorties  are  lost  forever  in  this 
case. 
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TABLE  4 


INPUT  VARIABLES  IN  TRAIL 


Name  of 
variable 

Symbol  for 
variable 
in  text 

t 

Definition 

DR 

r 

1 

i 

Reconnaissance  aircraft  radar  detection 
range 

DRH 

r 

Ttailerborne  helicopter  radar  detection  j. 

range 

i 

EPS 

i 

Tolerance  in  calculation  of  integral  in  1 

equation  28  (example:  .001  for  three- 

decimal-place  accuracy)  '■ 

FRE 

Frequency  vector  (weights  attached  to  each 
of  the  aircraft  mission  radii  in  R, 

according  to  the  likelihood  of  each  mission  , i 

radius.  Sums  to  1.) 

HTR 

To 

Trailerborne  Jielicopter  reaction  time —  1  1 
from  breaking  of  trail  to  initiation  of  | ; 
search 

N 

n 

jjfamber  of  days  trip  lasts 

NS 

*8 

1  i 

JAimber  of  jorties  a  reconnaissance  air¬ 
craft  makes  every  month  i 

NP 

ft’ 

|  ; 

Jjumber  of  ships  (to  be  trailed)  at  sea  1 

at  any  given  time  ^ 

OPP 

nopp 

Number  of  opportunities  to  break  trail  on  a 

day  when  circumstances  (e.g.,  weather)  are  | 

propitious  for  trail-breaking 

PC 

po 

(■ 

^Probability  ship  is  trailed  when  exiting  i 

port  ! 

PS 

p 

_Probability  of  finding  ship  by  recon-  j 

naissance  aircraft,  on  which  aircraft 
inventories  are  based 

R 

Vector  bearing  the  mission  radii  the  air¬ 
craft  might  be  called  on  to  cover 

i 
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TABLE  4  (Cont'd) 


Name  of 
variable 

Symbol  for 
variable 
in  text 

Definition 

S0 

so 

Number  of  hours  in  a  day  (»24) 

Si 

81 

Number  of  days  in  a  month  (=30) 

SMX 

^nax 

Time  from  trail-breaking  when  trailer- 
borne  helicopter  stops  searching  and 
calls  reconnaissance  aircraft 

TKF 

q 

Number  of  tankers  required  per  recon¬ 
naissance  aircraft 

TR 

- 

Time  from  trail-breaking  to  reconnais¬ 
sance  aircraft  taking  off  (TK>SMX) 

TS 

CS 

Reconnaissance  aircraft  time  on  station 
for  each  of  the  mission  radii  in  R 
(vector) 

V 

V 

Ship's  cruise  speed  (long-term  rate  of 
expansion  of  ship's  radius  of  uncertainty) 

VI 

V1 

Rate  of  expansion  of  ship's  radius  of 
uncertainty  in  directions  towards  the 
trailer,  shortly  after  breaking  trail 

V2 

v2 

Rate  of  expansion  of  ship's  radius  of 
uncertainty  in  directions  away  from  the 
trailer,  shortly  after  breaking  trail 
(usually,  ship's  maximum  sustained  speed) 

VA 

- 

Reconnaissance  aircraft  transit  speed  to 
station 

VH 

v' 

frailerborne  helicopter  search  speed 

VP 

v' 

Reconnaissance  aircraft  search  speed 

VT 

VT 

Jtailer-to-ship  closing  speed  after  air¬ 
craft  has  detected  ship 

WEA 

w 

Probability  that  circumstances  (e.g. 
weather)  will  be  propitious  for  trail¬ 
breaking  on  a  given  day . 
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TABLE  5 


OUTPUT  VARIABLES  IN  TRAIL 


Name  of 
variable 

Symbol  for 
variable 
in  text 

Definition 

Cl 

C1 

Quantity  defined  in  equation  9  (vector) 

C2 

c2 

Quantity  defined  in  equation  9 

E 

E 

Uncorrected  fraction  of  time  ship  is 
free  of  trail  (equation  29)  (vector) 

E0 

Eo 

Fraction  of  time  ship  is  free  of  trail 
when  P'  ■  0  (no  reacquisition  capability) 
(equation  39) 

EB 

E 

Corrected  fraction  of  time  ship  is  free 
of  trail  (equation  40a)  (vector) 

EBjO 

- 

Average  value  of  EB 

FBN 

f(n) 

Expected  number  of  trail-breaking  events 
during  ship's  n-day  trip  (equation  47) 
(vector) 

II 

- 

Value  of  integral  in  equation  28 
(vector) 

KK 

k 

Constant  such  that  trail  resumes  at  time 
kT  if  aircraft  detects  ship  at  time  T 
(equation  20) 

NR 

NR 

Required  ^lumber  of  reconnaissance  aircraft 
(equation  62)  (vector) 

NRT 

Nrt 

Total  number  of  aircraft  jrequired,  including 
tankers  (equation  64)  (vector) 

NRT0 

- 

Average  value  of  NRT 

NU 

V 

Number  of  trail-breaking  events  per  month 
that  is  not  exceeded  more  than  2  percent 
of  the  time  (equation  63)  (vector) 

TABLE  5  (Cont'd) 


UNCLASSIFIED 

Symbol  for 

Name  of  variable 

variable  in  text  _ Definition _ 

P  P  Effective  probability  ship  shakes  free  on  a 

day  when  it  starts  out  trailed  (equation  II) 
(vector) 

PH  -  probability  trallerborne  Jjelicopter 

reacquires  ship  per  trail -breaking 
attempt  (p0(Tmax)  in  equation  11) 

PMAX  Piaax  Maximum  probability  of  aircraft  detect¬ 

ing  ship  (given  unlimited  search) 

(equation  10)  (vector) 

PP  P*  Effective  probability  ship  is  acquired 

by  its  trailer  on  a  day  when  it  starts  out 
free  (equation  28)  (vector) 

QN  Qp  Probability  ship  is  free  on  day  n — the  last 

day  of  the  trip  (equation  36)  (vector) 

S  s  JSteady -state  daily  probability  of  either 

breaking  or  reacquiring  trail  (equation 
52)  (vector) 

T0  Tq  JElme  since  trail-breakirsg  when  recon¬ 

naissance  aircraft  arrives  on  station 
(vector) 

T1  tj  Aircraft  _time  on  station  per  trail-breaking 

event  that  guarantees  a  probability  PS  (p) 
of  finding  ship  («(KK  x  TAU)  -  T0)  (vector) 

TAU  t  Aircraft  search  time  that  results  in  a 

probability  PS  of  finding  the  ship 
(equation  24)  (vector) 

TT  tt  Reconnaissance  aircraft  _transit  _time  to 

station  (vector) 

AB,C11,E1,F,G, 

H,I,INT,K,L,Q,  Variables  used  internally 

W,X,X1,X2,Y 
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TABLE  6 


APL  CODING  OF  TRAIL 


9ARMCQ3* 

9  ARM  $  Q 
[13  *  INPUT*  • 

C23  0^*«-3+QIO«-l  #  Of  G«*0  •  98 

[33  '  VAsb  1  f  VA  J  •  p  TRs'jTRj',  VPe'fVPf*,  DRa  '  )  OR  }  1  r  Vs'fVf', 

psa*}ps;>,  vts'jvtj*,  so=,;so;,f  epsb'jeps;',  po=*}po 

i't  Na • fN| • ,  NPa'JNPJ*,  NSb'JNSJ*,  SI  a  *JS1J*,  Via** 
VU'r  V2=*»V2  »'»  VHs  •  *  VM  f  •  p  DRHa  '  f  DRH  f  '  f  HTRs  ’  }  HTR  f  '  f 
OPPa*|OPP}'f  SMXa •  | SMX |  *  f  ME As  *  j  WE A *  1  ,TKf  a  **TKF 

[43  '*  *  *>R 

C53  'FRE  a  ' f FRE 
C63  ’TS  a  ' ITSJ1  6/ 

[73  • output • | i  6f'-' 

[83  '<?2  =  •  ;c2<-(2xOAxVP)  j.(Vn2)xol 

[93  '««  ■  ' i KKFl+ViVTx2*0.5 

[103  'AVG  TT  a  1 }+/FREXTT; ' , TT  s  ' }TTfRiVA 
[113  'AVG  TO  *  1 ; +/FREXT0J * ,  TO  a  *JT0«-TT  +  TR 

[123  'Cl  a  •  ;ci«.«  (0-C2)^T0 

[133  1  AVG  PMAX  a  '  }+/FRExPMAX{  •  f  PMAX  a  «|PMAX«-1-C1 

C143  4<PS>L/FMAX)fX2 

[153  ‘AVG  TAU  a  1  }  +/FRE  X  TAU  J  '  ,  TAU  a  '  )  T  AUfC2r«  (  1  -PS  )  -j-C  1 

[163  IFIXFO 

C173  Xi:AB4.(iTAU[I3)r^.TO[I«.X  +  13 
C183  C11«-C1[I] 

[193  XX«-XX,AB  INTEGRAL  EPS 

C203  -»<i<pR>fXi 

[213  *AVG  II  a  • (+/FREXII J • ,  XX  a  *|Il4-l*II 

[223  'AVG  FP  ■  'J+/F RExPPJ',  PP  a  '  JPPI-1-C1XC2XXX  +  PS 

[233  fee 
[243  e i fP+pp 

[253  '«VG  E  a  1  t  +/FRE  XE  f  '  f  E  a  •  )  E «- (  P  + E 1  )  +  (  (  PP - PO  X  e  1  )  X  (  1 -C  1 
)X1-(1-E1)*N)^NXE1*2 

[263  'EO  a  '  ;eof1-fox(1~F)x(1-<1-p>*n)^.nxp 

[273  *«VG  EB  a  •  $ebo«-+/frexeb;  •  ,  EB  a  •  *EB«-<  (PSXE)+E0X1-PS) 
*G 

[283  ,ftvG  QN  =  1  I+/FREXON)  •  ,  ON  a  '  J  ON*.  (  P+  <  PP-PO  X  E  1  )  X  (  1  -E  1  ) 

*N)+el 

[293  'AVG  FBN  a  • f +/FRE x FBN } • f  FBN  a  • } FBNfPX ( NX 1-E > +P0-1 - 
ON 


-33- 


i iii  1 1  i'i  iiOiiiiiiiiMaiiMiiio 


TABLE  6  (ConC'd) 


C303  *AVO  *  ■  ,*+/AAEX*»,f  s  a  'fS«.AxAA*El 

C313  *avg  nu  *  •  }+/r*cxHUi  •  ,  hu  ■  •  »nu«-<s  j+n>  x  (naxabn)  +2*  <  I 
NAxNXSxl-S> *0.5 

C323  ‘AVO  T1  a  1  t  +/FREXT1  f  •  ,  TJ  a  *  >  T  1 «-  (  KK  X  TAU  ) -TO 
C333  'AVO  M*  *  *  )+/FRE!(HR)  •  f  HR  a  1  f  NRfNUX  T 1 -rNS  X  TS 

C343  'OVO  HUT  S  •  ;HRTO«-  +  /frKEXNRT|  •  ,  NAT  a  •  f  NAT*.  (  \  +T  KF  )  XMR  | 
1  1/'  * 

[353  'SUMMARY*  AVO  NAT  a  '  } NAT  0  f  '  »  AVO  EB  a  '{C»0;'r  *0  =  ' 

»«0»1  60t'-' 

C3  61  ■* 0 

C373  X2J  'CASE  INVALID,  AS  IS  UNATTAINABLE ,* %  1  60f ' ~  ' 

V 

VZNTEGAALC03? 

V  ZNT  f AB  ZNTEGAAL  EAS|N;W|K;L 

C13  **Alf0.5*(  <  M«-0)++/*HTEGRAND  A»)xHf-/|A» 

C23  *:««-«»  (N[  <f  N)  -  M-l  3-S-2)  +m*  +  /xnte°AAND(  1  fOB)  +  <H«.Ni"2)  *“1+2* 

\2R’*1«R4>K  +  X  +  L«.0 

C33  OJ«M-WrxNT<-(  (  (4*L)xNCfN3)-»»C{f»*)-it3)+-i  +  4*L«.L+i 
C43  -»OX  |  L  <  K 

C53  -»rx\EAS<  |-/wc<fM«.(-i4-2xK)tw)-0f  M3 

V 

VINTEGAAND[Q]V 

V  Yf INTEGRAND  X 

C13  T«-A((S0XC2^MM)XX*2)+XXC2+(S0+KI<)XBC11 

v 

BAECCQ3V 

V  AEEJG 

C13  o<-(4xDAHxvH)i0ix  <  vi*2)+v2*2 
C23  ’*«  *  '^Htl-nOKiSMXJ-iHTR 

C33  *  '  |a«-weax1-ah*oaa 

V 
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integral  be  too  small,  space  limits  may  be  reached,  and  EPS  may  have  to 
be  reduced  (which  sometimes  happens  for  small  ship  speeds  v  because  of 
the  large  value  of  the  integral). 


Several  aircraft  mission  radii  can  be  input  at  the  same  time,  with 
their  relative  frequencies  of  occurrence.  The  results  are  then  pre¬ 
sented  separately  per  mission  radius  and  averaged. 

A  word  on  units.  The  input  and  output  lists  refer  to  hours,  days, 
and  months.  This  choice  is  arbitrary  and  can  be  changed.  The  issue 
here  is  not  preference  but  reasonableness  of  results.  A  ship  is  in  the 
“trail"  or  "no-trail"  state  for  a  whole  day.  If  it  is  felt,  for  exam¬ 
ple,  that  a  12-hour  interval  would  be  more  appropriate  than  a  day;  the 
input  S0  must  be  changed  from  24  to  12.  The  same  goes  for  the  month, 
which  is  the  time  interval  on  which  aircraft  Inventories  are  based 
(input  SI). 

Examples 

A  sample  run  of  TRAIL  is  shown  in  table  7.  The  input  quantities 
are  arbitrary  and  are  not  Intended  to  reflect  specific  systems. 

An  example  of  the  kind  of  results  TRAIL  can  generate  is  the  frac¬ 
tion  of  ships  at  sea  that  are  free  of  trail  at  any  given  time,  E,  as  a 
function  of  the  number  of  aircraft  the  enemy  deploys.  This  function  is 
obtained  by  running  TRAIL  several  times  for  different  values  of  p  (the 
trail  reacquisition  probability  that  the  enemy  tailors  his  aircraft 


inventory 

to  achieve). 

In  our 

illustrative 

case. 

the 

result 

is: 

P 

(PS) 

0 

.1 

.3 

.5 

.7 

.8 

.9 

.95 

.99 

S" 

(NRT) 

0 

12 

18 

22 

26 

29 

36 

45 

106 

(EB) 

.84 

.82 

.68 

.52 

.35 

.27 

.18 

.14 

.11 

Figure  9  is  a  plot  of  E  versus  nR,j,.  When  no  aircraft  are  deployed, 
the  trailing  ships  are  left  to  their  own  devices. 


TABLE  7 


SAMPLE  RUN  OF  TRAIL 


ap-M 
T  HPUT  • 

tf-5,  vp-150  ,  **=250,  v=r20,  *=-s  =  o.7>  vt  =  3,  s<>=24,  eps 
-  1  e  “  5 ,  p'0=0  •  5 »  M= 120,  ”*=50,  "S^IO,  si  «  30,  vi  =  io,  v 

2-20»  vM,=  75f  orm=4o,  htr  =  i.2,  npp=4,  sm><  =  3»  weo=o.05  »t 

kp  s  0.9 

*  *  2000  3000  4000 

*■**  -  0.2  0.7  0.1 

T"  25  15  10 


n  U  T  P  l' t 


-'2  „  59.69 
**  '  5.714 

«'•’<*  TT  -r  4 . 462 » TT  -  3.077  4.615  6.154 

ovr,  T0  =  9 . 462  »  t  0  =  8.077  9.615  11.15 

ri  -  0 . 0006178  0.002015  0.004744 

rtvr,  PMfl =  0.998,  =  0.9994  0.998  0.9953 

av<i  Tfiu  r:  H.72,  Tflu  =  9. £49  H.93  14.39 

«ve  it  ■  6.674,  ii  =  15.29  4.862  2.122 

avo  pr  =  0.1683,  ff  =  0.1946  0.1646  0.1417 

fh  »  0.8519 

f  -  0.02367 

ovr.  e  »  0.1372,  e  =  0.1201  0.1392  0.1581 

r 0  =  0.8378 


AVG 

EG 

*  0.3545, 

EB 

=  0.3423  0.3559  0.3694 

AVG 

OH 

=  0.124, 

OH  = 

0.1084  0.1257  0.1431 

AVG 

FBI) 

-  2.441, 

FBM 

*  2.49  2.436  2.383 

AVG 

S  sr 

0.02073, 

S  -- 

0.0211  0.02069  0.02028 

AVG 

MU 

*  36.03, 

MU  s 

36.69  35.96  35.24 

AVG 

T1 

*  57.5,  ti  « 

47.06  58.55  71.08 

A  V  G 

MR 

~  13.71, 

MR  •= 

6.905  14.04  25.05 

AVG 

MRT 

=  26.05, 

MRT 

=  13.12  26.67  47.6 

summary «  AVG  M ft t  s  26.05,  «VG  EB  =  0.3545,  EO  =  0.8378 
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Fraction  of  fleet  free  of  trail  IE) 


FIG.  9:  SAMPLE  RESULTS  OBTAINED  FROM  TRAIL  RUNS8 


* Baud  on  a  fleet  of  50  target  ships  and  50  trailing  ships  at  sea  all  the  time. 
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